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Parvan [arXiv:0911.0383V l] [3| has recently presented some calculations in order to demonstrate 
the incorrectness of the results obtained from the generalized multinomial coefficients (GMC) pre- 
sented in Ref. Q]. According to Parvan, the aforementioned approach of studying maximum entropy 
probability distributions is erroneous. In this comment 1 demonstrate that Parvan's arguments do 
not hold true and that the obtained results from GMC do not present either mathematical or 
physical discrepancies. 

D ' PACS numbers: PACS: 05.: 05.20.-y: 05.30.-d 

' The fundamental flaw Parvan believes to have discovered in Ref. is exactly the same misunderstanding, with 
which E. T. Jaynes had once been acquainted The subject matter is the difference between the Boltzmann and 
I— !• Gibbs entropies. In order to shed light on this issue, let us define W as the number of the allowed energy states that a 
\ particle can occupy within an ensemble, and pi are the respective probabilities (or simply called particle probabilities) 
^ ■ for each one of these particles. Then, the Boltzmann entropy reads 

S' 

> 

I ; , where N is the size of the collection (number of particles). Denoting the total number of configurations associated 

d • with the energy states by fl (F in the notation of Ref. [i|), we then assign the occurrence probability Pj to each 

^ [ configuration. This results in the following expression for the Gibbs entropy Sq 

g: SG = Y.p,Hyp,)- (2) 

,o,. 

with f2 = {Iq. The two expressions given above are not always equal to each other Q. In fact, their equality is 
CS| ■ satisfied only when the particles are independent from one another 0]. If this is the case, then the total number 
^ ' of configurations in the Gibbs ensemble equals to that of the Boltzmann entropy i.e., ilo = ^^bg- This property 
25 ' uniquely characterizes the BG-ensemble, and clearly distinguishes it from the Bose-Einstein-Gibbs and Fermi-Dirac- 
Gibbs ensembles. Accordingly, one can use Sb or Sq interchangeably in BG-collection of configurations. From the 
aforementioned equality, one can directly verify the extensivity of the BG-entropy. It is also well known, that from 
the combinatorial point of view Sb can be derived on the basis of the multinomial coefficient 
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Sb = NJ2P'M^/P'd, (1) 



in the limit — > cxd. Then, for equal probabilities one obtains the BG-configuration function depending on N, namely 

U ■ ^BoiN) = lim C^^ = exp[Nln{W)) . (4) 

In order to find the thermodynamic expression of ^Ibg however, one should replace W with physical quantities, since 
the latter is of purely combinatorial nature. 

Taking above explanation into consideration, one may think to generalize Eq. (|3]) in order to derive a generalized 
entropic form. Indeed, Suyari first studied such a possibility in Ref. Q with Tsallis (T) entropy being the point 
of departure. His success was partial in the sense that one could clearly see the connection between his generalized 
multinomial coefficient (GMG) and T-entropy but its exact structure could not be derived for all values of the 
generalization parameter. Following Suyari's steps, I continued along the same direction but in a different way in 
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Ref. so as to define two CMC's (CMCi and CMC2) in order to derive (among others) (Boltzmann-)Tsallis and 
(Boltzmann-)Renyi entropic structure and study their maximum probability distributions. Since the aforementioned 
entropic structures concern the generahzation of the BC-ensemblc, whether ones uses Boltzmann generalization or 
Cibbs generalization does not matter. Even from this much explanation, it becomes evident that Parvan's doubt of 
considering the entropy definitions in Ref. Q as incorrect is misplaced. 

Let us now recapitulate the results obtained in Ref. Q . Regarding Tsallis entropy, it was shown that the MaxEnt 
distribution structure obtained from CMCi is in accordance with the one obtained from Jaynes Formalism (JF) 
with ordinary linear constraints, 1/ expg(|x|) with q G [0, 1], while the MaxEnt distribution structure obtained from 
CMC2 is in accordance with the ones obtained from JF with escort linear constraints, 1/ exp2_^(|x|) with q G [1, 00). 
The g-intervals stem from the requirement of the positivity of the CMC's, since a negative number of states is 
senseless. Unfortunately, reading again Ref. [2] after Parvan's critique, I realized that, although the equations 
themselves demonstrate that in Tsallis case (and only in this case) the results from JF and CMC's coincide, in 
the text right after Eq. (61) in Ref. [2], it states that all the results obtained from CMC's are in contradiction 
with the ones obtained from JF. This is a misprint, which might be the source of confusion on part of Parvan. 
This sentence only refers to Rcnyi and Caussian entropy as it becomes obvious in the context later, where it is 
demonstrated that the combination of the results obtained from both methods. JF and CMC, sheds light on the vi- 
olation of concavity of Tsallis entropy in its Escort Distribution Representation for q e [0, 1], also discussed in Ref. [H. 

According to the results in Ref. (H the configuration function, or following Parvan's notation, the total number of 
microstatcs in Tsallis and Renyi ensemble are of the following form 

flrr = expq(x), JIr = exp(x) , (5) 

which are in agreement with the respective entropy extensivity 0. Parvan, arbitrarily and erroneously, chooses 
to adopt Eqs. (27)-(29) in his critique i.e., the same V, for all three ensembles, BC-, T- and R-ensemble. It is 
then straightforward that replacing for example a configuration function of ordinary exponential type in the Tsallis 
entropy for equal probabilities, one violates the extensivity property, as in Eq. (28) in Ref. [![. 

Regarding Rcnyi entropy, it is far more obvious that since the averaging procedure in Rcnyi entropy is of exponential 
type and not of linear type as in BC-entropy i.e., 

^ = (-^)nn = f:P^n, I = (-.)e.p = In (f:P^e^'-^A (6) 

i=l \i=l / 

where ln(l/pi), the entire Renyi ensemble is characterized by the aforementioned non-linear averaging procedure. 
In fact, historically, Renyi obtained his measure by conforming to the requirement of additivity. The preservation of 
additivity leaves one with only two choices if one would like to keep the definition of information gain undeformed. 
These two choices are linear averaging and exponential averaging procedures. The former paves the way to BC 
entropy, whereas the latter yields to Rcnyi measure. In this sense, what Rcnyi did is remarkable, since one cannot 
write a new undeformed generalized entropy measure. Renyi saw the second possibility and that is it. However, 
it was not his intention to write a generalized entropy measure valid for a generalized Thermostatistics. This is a 
different agenda then just generalizing BC measure. Thus, in a thermostatistical framework, whenever one applies 
JF, or equivalently Parvan's formalism Q on Renyi entropy, one should consider the proper constraints, as correctly 
noticed in Ref. (8|. Then, the derived MaxEnt probability distribution is of ordinary exponential type and it results 
simultaneously and consistently the correct interval of concavity for the deformation parameter associated with the 
Sr It is worth noticing that the same g-range is obtained through CMC-method. It is indeed surprising that the 
maximization of Renyi entropy so far was carried out through linear averaged constraints. However, one can even 
obtain inverse power law stationary distributions from the maximization of BC entropy with arbitrary constraints if 
one is willing to do arbitrary physics. For serious physicists though, this is not an option. 

Further, I would like to present some more evidence about the correctness of CMC-approach following a different 
mathematical route. The standard method in computing MaxEnt distributions is by applying JF to a functional 
including the entropy measure and the suitably averaged constraints. In case of BC-statistics, the results obtained 
from the former formalism coincide with the ones derived within thermostatistical theory. However, once we apply JF 
on generalized entropic structures, there are no reference results to compare with and assure their correctness. Indeed, 
the probabilities are mathematically dependent quantities through the normalization condition, thus the introduction 
of this dependence as a constant in JF assuming initially mathematical independence between the probabilities may 
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lead generally to wrong results. So, the idea is to check which probability distribution structure for a given entropy 
expression is compatible with the fundamental equilibrium relation of Thermodynamics, namely 

dQ = dE + dW, (7) 

where Q, E and W is the heat flow, the internal energy and the produced work of the physical system under 
consideration, without computing derivatives. If we can determine the aforementioned structure, we can conclude 
that even the MaxEnt distribution must be of the same form. If we choose the initial values equal to zero, the above 
equation is written as 

Q = E + W. (8) 
The statistical representation of this equation of ordinary statistics is given as 

^MV^O)^, = (^.>ii„ + (W)ii„ =^ ^iln(l/P,)-e,->V^^ =0, (9) 

From the right hand relation of Eq. ([S]), it becomes obvious that in the ordinary case the probabilities which are 
compatible with the BG-cntropic structure and the thermodynamic relation ([7)) are of the following structure 

l/P, = hi-i (/3e, + /3 W) = cxp (p (e, + W) ) . (10) 

This is the probability distribution obtained within the ordinary theory of Thermodynamics. Indeed, if we replace 

the above probability distribution in the probabilistic BG-heat expression Q = ( -g ln(l/-Pi)) , we obtain again the 

I— — » 

fundamental equilibrium thermodynamic relation in Eq. ([8|). Comparing Eqs. (|10p and (|4]) we see that the statistical 
quantity W takes the form W = exp (-^(ej + W)) in terms of physical quantities. 

Considering the heat expression Q, we observe two basic possibilities to generalize it, either we generalize the 
uncertainty ln(l/Pi) — > ln^(l/i-i) or the averaging procedure (• • •)[jjj — >■ (• • •)^ (or both of them), where ^ = {Ci}i=i,...,ti 
and a = {(Ti}i=i,...,„ are two deformation parameter sets. For certain values of these parameters ^ — and a — >■ tro 
the generalized expression recovers the ordinary one in Eq. ([9]) . In the first case of generalized uncertainty (trace- form 
entropies) Eq. (jH]) tends to 

1 . A /I 



^Inai/P,)^^ = (e»)iin + (W)ii„ =^ ^-lne(l/P,)-£.-W^^ =0 (11) 

and the obtained probability distribution is given by 

l/P, = cxp^ (/3 (e, + W) ) . (12) 

The distribution structure in Eq. ((T2|) is the one obtained from GMCi in Rcf. ^2] for Tsallis entropy. We also see 
here that all trace-form entropies are compatible with linear averaged physical quantities. Pi in Eq. (|12p preserves 
the extensivity property of the respective entropy, as explained in Ref. [6| . 

In the second case of the generalized nonlinear average (non-trace-form entropies) Eq. (|9]) tends to 

1 



^ln(l/PO^ + (13) 

with (• • •)g._>cro ~ (' ■ ')iin- '^^^ determination of the compatible probability distribution cannot be done in the general 

case, thus each non-trace form entropy should be studied separately. For the Rcnyi entropy, a = exp, Eq. ((T^ takes 
the form 

{-pHllPr)) =(e.)exp + (WUp =^ ^P,^''"''-exp((l-g)(e, + W))^^ =0. (14) 

Then, the compatible probability distribution is of the form 

1/P, = exp (/? [e, + W) ) . (15) 
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As can be seen, this is exactly the distribution one obtains with the GMC-approach in Ref. [2]. Herewith, we 
complete the proof of the correctness of the results presented in Ref. , which are also confirmed by the general and 
consistent (contrary to the statement by Parvan who considers this maximization as a special case) maximization of 
Renyi entropy given in Ref. [1]. A crucial issue in the Eqs. (P))- ([T5|) is whether the quantities Si and W are the same 
in all ensembles. However, the answer to this question is not within the scope of this comment and deserves to be 
examined elsewhere in detail. 

The author would like to thank U. Tirnakli for bringing Ref. [l| to my attention. G. B. Bagci and U. Tirnakli are 
acknowledged for carefully reading this comment and fruitful discussions. 
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